Trends, patterns and relationship of antimicrobial use and resistance in bacterial isolates tested between 2015-2020 in a national referral hospital of Zambia.

Increased antimicrobial resistance (AMR) among bacteria underscores the need to strengthen AMR surveillance and promote data-based prescribing. To evaluate trends and associations between antimicrobial usage (AMU) and AMR, we explored a dataset of 34,672 bacterial isolates collected between 2015 and 2020 from clinical samples at the University Teaching Hospital (UTH) in Lusaka, Zambia. The most frequently isolated species were Escherichia coli (4,986/34,672; 14.4%), Staphylococcus aureus (3,941/34,672; 11.4%), and Klebsiella pneumoniae (3,796/34,672; 10.9%). Of the 16 drugs (eight classes) tested, only amikacin and imipenem showed good (> 50%) antimicrobial activity against both E. coli and K. pneumoniae, while nitrofurantoin was effective only in E. coli. Furthermore, 38.8% (1,934/4,980) of E. coli and 52.4% (2,079/3,791) of K. pneumoniae isolates displayed multidrug resistance (MDR) patterns on antimicrobial susceptibility tests. Among S. aureus isolates, 44.6% (973/2,181) were classified as methicillin-resistant (MRSA). Notably, all the MRSA exhibited MDR patterns. The annual hospital AMR rates varied over time, while there was a weak positive relationship (r = 0.38, 95% CI = 0.11-0.60) between the monthly use of third-generation cephalosporins (3GCs) and 3GC resistance among Enterobacterales. Overall, the results revealed high AMR rates that fluctuated over time, with a weak positive relationship between 3GC use and resistance. To our knowledge, this is the first report to evaluate the association between AMU and AMR in Zambia. Our results highlight the need to strengthen antimicrobial stewardship programs and optimize AMU in hospital settings.

Detection of DNA Modification Using Nanopore Sequencers.

DNA modification is a crucial factor of epigenetic modification and has vital functions for gene regulation and phenotype control. A profound understanding of DNA modification requires precise mapping of the modified bases on genomic DNA. In addition to methods such as bisulfite sequencing and single-molecule real-time (SMRT) sequencing of PacBio sequencers, nanopore sequencers can be also utilized for the detection of DNA modification. Here, I will briefly review the three methods for the detection of DNA modification with nanopore sequencers and introduce a protocol using MinION and Megalodon.

Identification, expression, and purification of DNA cytosine 5-methyltransferases with short recognition sequences.

BACKGROUND: DNA methyltransferases (MTases) are enzymes that induce methylation, one of the representative epigenetic modifications of DNA, and are also useful tools for analyzing epigenomes. However, regarding DNA cytosine 5-methylation, MTases identified so far have drawbacks in that their recognition sequences overlap with those for intrinsic DNA methylation in mammalian cells and/or that the recognition sequence is too long for fine epigenetic mapping. To identify MTases with short recognition sequences that never overlap with the CG dinucleotide, we systematically investigated the 25 candidate enzymes identified using a database search, which showed high similarity to known cytosine 5-MTases recognizing short sequences. RESULTS: We identified MTases with six new recognition sequences, including TCTG, CC, CNG, TCG, GCY, and GGCA. Because the recognition sequence never overlapped with the CG dinucleotide, MTases recognizing the CC dinucleotide were promising. CONCLUSIONS: In the current study, we established a procedure for producing active CC-methylating MTases and applied it to nucleosome occupancy and methylome sequencing to prove the usefulness of the enzyme for fine epigenetic mapping. MTases that never overlap with CG dinucleotides would allow us to profile multiple epigenomes simultaneously.

Chlamydia trachomatis L2/434/Bu Favors Hypoxia for its Growth in Human Lymphoid Jurkat Cells While Maintaining Production of Proinflammatory Cytokines.

The role of lymphocytes as a cornerstone of the inflammatory response in the invasive pathogenesis of Chlamydia trachomatis (Ct) LGV (L1-3) infection is unclear. Therefore, we assessed whether the adaptation of CtL2 to immortal lymphoid Jurkat cells under hypoxic conditions occurred through proinflammatory cytokine profile modification. The quantities of inclusion-forming units with chlamydial 16S rDNA confirmed that CtL2 grew well under hypoxic rather than normoxic conditions in the cells. Confocal microscopic imaging and transmission electron microscopy revealed the presence of bacterial progeny in the inclusions and showed that the inclusions were larger under hypoxic rather than normoxic conditions; this was supported by the results of 3D image construction. Furthermore, PCR-based analysis of proinflammatory cytokines revealed that the gene expression levels under hypoxic conditions were significantly higher than those under normoxic conditions. In particular, the expression of two genes (CXCL8 and CXCR3) was significantly diminished under normoxic conditions. Taken together, the results indicated that hypoxia promoted CtL2 growth in Jurkat cells while maintaining the levels of proinflammatory cytokines. Thus, Ct LGV infection in lymphocytes under hypoxic conditions might be crucial to a complete understanding of the invasive pathogenesis.

Chlamydia trachomatis Requires Functional Host-Cell Mitochondria and NADPH Oxidase 4/p38MAPK Signaling for Growth in Normoxia.

Chlamydia trachomatis (Ct) is an intracellular energy-parasitic bacterium that requires ATP derived from infected cells for its growth. Meanwhile, depending on the O2 concentration, the host cells change their mode of ATP production between oxidative phosphorylation in mitochondria (Mt) and glycolysis; this change depends on signaling via reactive oxygen species (ROS) produced by NADPH oxidases (NOXs) as well as Mt. It has been proposed that Ct correspondingly switches its source of acquisition of ATP between host-cell Mt and glycolysis, but this has not been verified experimentally. In the present study, we assessed the roles of host-cell NOXs and Mt in the intracellular growth of CtL2 (L2 434/Bu) under normoxia (21% O2) and hypoxia (2% O2) by using several inhibitors of NOXs (or the downstream molecule) and Mt-dysfunctional (Mtd) HEp-2 cells. Under normoxia, diphenyleneiodonium, an inhibitor of ROS diffusion, abolished the growth of CtL2 and other Chlamydiae (CtD and C. pneumoniae). Both ML171 (a pan-NOX inhibitor) and GLX351322 (a NOX4-specific inhibitor) impaired the growth of CtL2 under normoxia, but not hypoxia. NOX4-knockdown cells diminished the bacterial growth. SB203580, an inhibitor of the NOX4-downstream molecule p38MAPK, also inhibited the growth of CtL2 under normoxia but not hypoxia. Furthermore, CtL2 failed to grow in Mtd cells under normoxia, but no effect was observed under hypoxia. We conclude that under normoxia, Ct requires functional Mt in its host cells as an ATP source, and that this process requires NOX4/p38MAPK signaling in the host cells. In contrast to hypoxia, crosstalk between NOX4 and Mt via p38MAPK may be crucial for the growth of Ct under normoxia.

Cas9-guided haplotyping of three truncation variants in autosomal recessive disease.

An autosomal recessive disease is caused by biallelic loss-of-function mutations. However, when more than two disease-causing variants are found in a patient's gene, it is challenging to determine which two of the variants are responsible for the disease phenotype. Here, to decipher the pathogenic variants by precise haplotyping, we applied nanopore Cas9-targeted sequencing (nCATS) to three truncation COL7A1 variants detected in a patient with recessive dystrophic epidermolysis bullosa (EB). The distance between the most 5' and 3' variants was approximately 19 kb at the level of genomic DNA. nCATS successfully demonstrated that the most 5' and 3' variants were located in one allele while the variant in between was located in the other allele. Interestingly, the proband's mother, who was phenotypically intact, was heterozygous for the allele that harbored the two truncation variants, which could otherwise be misinterpreted as those of typical recessive dystrophic EB. Our study highlights the usefulness of nCATS as a tool to determine haplotypes of complicated genetic cases. Haplotyping of multiple variants in a gene can determine which variant should be therapeutically targeted when nucleotide-specific gene therapy is applied.

A GCDGC-specific DNA (cytosine-5) methyltransferase that methylates the GCWGC sequence on both strands and the GCSGC sequence on one strand.

5-Methylcytosine is one of the major epigenetic marks of DNA in living organisms. Some bacterial species possess DNA methyltransferases that modify cytosines on both strands to produce fully-methylated sites or on either strand to produce hemi-methylated sites. In this study, we characterized a DNA methyltransferase that produces two sequences with different methylation patterns: one methylated on both strands and another on one strand. M.BatI is the orphan DNA methyltransferase of Bacillus anthracis coded in one of the prophages on the chromosome. Analysis of M.BatI modified DNA by bisulfite sequencing revealed that the enzyme methylates the first cytosine in sequences of 5'-GCAGC-3', 5'-GCTGC-3', and 5'-GCGGC-3', but not of 5'-GCCGC-3'. This resulted in the production of fully-methylated 5'-GCWGC-3' and hemi-methylated 5'-GCSGC-3'. M.BatI also showed toxicity when expressed in E. coli, which was caused by a mechanism other than DNA modification activity. Homologs of M.BatI were found in other Bacillus species on different prophage like regions, suggesting the spread of the gene by several different phages. The discovery of the DNA methyltransferase with unique modification target specificity suggested unrevealed diversity of target sequences of bacterial cytosine DNA methyltransferase.

Complete Genome Sequences of 24 Strains of Bacillus cereus Isolated from Nosocomial Infection and Bacteremia Cases in Japan.

Bacillus cereus is mainly associated with foodborne illness but sometimes causes nosocomial infections. We previously reported that B. cereus strains of a specific sequence type, ST1420, were associated with nosocomial infection. Here, we determined the complete genome sequences of B. cereus strains isolated from nosocomial infection cases in Japanese hospitals.

Increased energy demand from anabolic-catabolic processes drives ß-lactam antibiotic lethality.

ß-Lactam antibiotics disrupt the assembly of peptidoglycan (PG) within the bacterial cell wall by inhibiting the enzymatic activity of penicillin-binding proteins (PBPs). It was recently shown that ß-lactam treatment initializes a futile cycle of PG synthesis and degradation, highlighting major gaps in our understanding of the lethal effects of PBP inhibition by ß-lactam antibiotics. Here, we assess the downstream metabolic consequences of treatment of Escherichia coli with the ß-lactam mecillinam and show that lethality from PBP2 inhibition is a specific consequence of toxic metabolic shifts induced by energy demand from multiple catabolic and anabolic processes, including accelerated protein synthesis downstream of PG futile cycling. Resource allocation into these processes is coincident with alterations in ATP synthesis and utilization, as well as a broadly dysregulated cellular redox environment. These results indicate that the disruption of normal anabolic-catabolic homeostasis by PBP inhibition is an essential factor for ß-lactam antibiotic lethality.

Detection of revertant mosaicism in epidermolysis bullosa through Cas9-targeted long-read sequencing.

Revertant mosaicism (RM) is a phenomenon in which inherited mutations are spontaneously corrected in somatic cells. RM occurs in some congenital skin diseases, but genetic validation of RM in clinically revertant skin has been challenging, especially when homologous recombination (HR) is responsible for RM. Here, we introduce nanopore Cas9-targeted sequencing (nCATS) for identifying HR in clinically revertant skin. We took advantage of compound heterozygous COL7A1 mutations in a patient with recessive dystrophic epidermolysis bullosa who showed revertant skin spots. Cas9-mediated enrichment of genomic DNA (gDNA) covering the two mutation sites (>8 kb) in COL7A1 and subsequent MinION sequencing successfully detected intragenic crossover in the epidermis of the clinically revertant skin. This method enables the discernment of haplotypes of up to a few tens of kilobases of gDNA. Moreover, it is devoid of polymerase chain reaction amplification, which can technically induce recombination. We, therefore, propose that nCATS is a powerful tool for understanding complicated gene modifications, including RM.

Clonal relationship between multidrug-resistant Escherichia coli ST69 from poultry and humans in Lusaka, Zambia.

Multidrug-resistant (MDR) Escherichia coli in food animals such as chickens is an emerging public health concern in Zambia. Additionally, the country's high demand for poultry products necessitates further investigation into the link between poultry and human MDR E. coli. Twenty cefotaxime-resistant E. coli isolates collected from poultry in Lusaka, Zambia, were screened for multidrug resistance and sequenced on MiSeq and MinION platforms. Genomes were assembled de novo and compared with 36 previously reported cefotaxime-resistant E. coli isolates from inpatients at the University Teaching Hospital, Lusaka. All (20/20, 100%) poultry isolates exhibited resistance to ampicillin, chloramphenicol and doxycycline. Phylogenetic analysis and hierarchical clustering showed a high degree of genetic relatedness between E. coli O17:H18-ST69 from poultry and humans. The E. coli O17:H18-ST69 clone accounted for 4/20 (20%) poultry- and 9/36 (25%) human-associated isolates that shared two plasmids harboring 14 antimicrobial resistance (AMR) genes. However, comparison analysis showed that the isolates also had other AMR plasmids distinct for each niche. Our results suggested clonal transmission of MDR E. coli between poultry and humans, with the potential acquisition of niche-specific AMR plasmids. Thus, the control of MDR E. coli requires a One Health approach involving both human and animal health sectors.

Risk factors and spatio-temporal patterns of livestock anthrax in Khuvsgul Province, Mongolia.

Anthrax is a worldwide zoonotic disease. Anthrax has long been a public health and socio-economic issue in Mongolia. Presently, there is no spatial information on carcass burial sites as a potential hazard of future anthrax outbreaks and possible risk factors associated with anthrax occurrences in Mongolia. Here, we analyze retrospective data (1986-2015) on the disposal sites of livestock carcasses to describe historical spatio-temporal patterns of livestock anthrax in Khuvsgul Province, which showed the highest anthrax incidence rate in Mongolia. From the results of spatial mean and standard deviational ellipse analyses, we found that the anthrax spatial distribution in livestock did not change over the study period, indicating a localized source of exposure. The multi-distance spatial cluster analysis showed that carcass sites distributed in the study area are clustered. Using kernel density estimation analysis on carcass sites, we identified two anthrax hotspots in low-lying areas around the south and north regions. Notably, this study disclosed a new hotspot in the northern part that emerged in the last decade of the 30-year study period. The highest proportion of cases was recorded in cattle, whose prevalence per area was highest in six districts (i.e., Murun, Chandmani-Undur, Khatgal, Ikh-Uul, Tosontsengel, and Tsagaan-Uul), suggesting that vaccination should prioritize cattle in these districts. Furthermore, size of outbreaks was influenced by the annual summer mean air temperature of Khuvsgul Province, probably by affecting the permafrost freeze-thawing activity.

Development of ELISA based on Bacillus anthracis capsule biosynthesis protein CapA for naturally acquired antibodies against anthrax.

Anthrax is a zoonotic disease caused by the gram-positive spore-forming bacterium Bacillus anthracis. Detecting naturally acquired antibodies against anthrax sublethal exposure in animals is essential for anthrax surveillance and effective control measures. Serological assays based on protective antigen (PA) of B. anthracis are mainly used for anthrax surveillance and vaccine evaluation. Although the assay is reliable, it is challenging to distinguish the naturally acquired antibodies from vaccine-induced immunity in animals because PA is cross-reactive to both antibodies. Although additional data on the vaccination history of animals could bypass this problem, such data are not readily accessible in many cases. In this study, we established a new enzyme-linked immunosorbent assay (ELISA) specific to antibodies against capsule biosynthesis protein CapA antigen of B. anthracis, which is non-cross-reactive to vaccine-induced antibodies in horses. Using in silico analyses, we screened coding sequences encoded on pXO2 plasmid, which is absent in the veterinary vaccine strain Sterne 34F2 but present in virulent strains of B. anthracis. Among the 8 selected antigen candidates, capsule biosynthesis protein CapA (GBAA_RS28240) and peptide ABC transporter substrate-binding protein (GBAA_RS28340) were detected by antibodies in infected horse sera. Of these, CapA has not yet been identified as immunoreactive in other studies to the best of our knowledge. Considering the protein solubility and specificity of B. anthracis, we prepared the C-terminus region of CapA, named CapA322, and developed CapA322-ELISA based on a horse model. Comparative analysis of the CapA322-ELISA and PAD1-ELISA (ELISA uses domain one of the PA) showed that CapA322-ELISA could detect anti-CapA antibodies in sera from infected horses but was non-reactive to sera from vaccinated horses. The CapA322-ELISA could contribute to the anthrax surveillance in endemic areas, and two immunoreactive proteins identified in this study could be additives to the improvement of current or future vaccine development.

Silkworm model for Bacillus anthracis infection and virulence determination.

Bacillus anthracis is an obligate pathogen and a causative agent of anthrax. Its major virulence factors are plasmid-coded; however, recent studies have revealed chromosome-encoded virulence factors, indicating that the current understanding of its virulence mechanism is elusive and needs further investigation. In this study, we established a silkworm (Bombyx mori) infection model of B. anthracis. We showed that silkworms were killed by B. anthracis Sterne and cured of the infection when administered with antibiotics. We quantitatively determined the lethal dose of the bacteria that kills 50% larvae and effective doses of antibiotics that cure 50% infected larvae. Furthermore, we demonstrated that B. anthracis mutants with disruption in virulence genes such as pagA, lef, and atxA had attenuated silkworm-killing ability and reduced colonization in silkworm hemolymph. The silkworm infection model established in this study can be utilized in large-scale infection experiments to identify novel virulence determinants and develop novel therapeutic options against B. anthracis infections.

Direct Regulons of AtxA, the Master Virulence Regulator of Bacillus anthracis.

AtxA, the master virulence regulator of Bacillus anthracis, regulates the expression of three toxins and genes for capsule formation that are required for the pathogenicity of B. anthracis. Recent transcriptome analyses showed that AtxA affects a large number of genes on the chromosome and plasmids, suggesting a role as a global regulator. However, information on genes directly regulated by AtxA is scarce. In this work, we conducted genome-wide analyses and cataloged the binding sites of AtxA in vivo and transcription start sites on the B. anthracis genome. By integrating these results, we detected eight genes as direct regulons of AtxA. These consisted of five protein-coding genes, including two of the three toxin genes, and three genes encoding the small RNAs XrrA and XrrB and a newly discovered 95-nucleotide small RNA, XrrC. Transcriptomes from single-knockout mutants of these small RNAs revealed changes in the transcription levels of genes related to the aerobic electron transport chain, heme biosynthesis, and amino acid metabolism, suggesting their function for the control of cell physiology. These results reveal the first layer of the gene regulatory network for the pathogenicity of B. anthracis and provide a data set for the further study of the genomics and genetics of B. anthracis. IMPORTANCE Bacillus anthracis is the Gram-positive bacterial species that causes anthrax. Anthrax is still prevalent in countries mainly in Asia and Africa, where it causes economic damage and remains a public health issue. The mechanism of pathogenicity is mainly explained by the three toxin proteins expressed from the pXO1 plasmid and by proteins involved in capsule formation expressed from the pXO2 plasmid. AtxA is a protein expressed from the pXO1 plasmid that is known to upregulate genes involved in toxin production and capsule formation and is thus considered the master virulence regulator of B. anthracis. Therefore, understanding the detailed mechanism of gene regulation is important for the control of anthrax. The significance of this work lies in the identification of genes that are directly regulated by AtxA via genome-wide analyses. The results reveal the first layer of the gene regulatory network for the pathogenicity of B. anthracis and provide useful resources for a further understanding of B. anthracis.

Novel chromosomal insertions of ISEcp1-blaCTX-M-15 and diverse antimicrobial resistance genes in Zambian clinical isolates of Enterobacter cloacae and Escherichia coli.

BACKGROUND: The epidemiology of extended-spectrum ß-lactamases (ESBLs) has undergone dramatic changes, with CTX-M-type enzymes prevailing over other types. blaCTX-M genes, encoding CTX-M-type ESBLs, are usually found on plasmids, but chromosomal location is becoming common. Given that blaCTX-M-harboring strains often exhibit multidrug resistance (MDR), it is important to investigate the association between chromosomally integrated blaCTX-M and the presence of additional antimicrobial resistance (AMR) genes, and to identify other relevant genetic elements. METHODS: A total of 46 clinical isolates of cefotaxime-resistant Enterobacteriaceae (1 Enterobacter cloacae, 9 Klebsiella pneumoniae, and 36 Escherichia coli) from Zambia were subjected to whole-genome sequencing (WGS) using MiSeq and MinION. By reconstructing nearly complete genomes, blaCTX-M genes were categorized as either chromosomal or plasmid-borne. RESULTS: WGS-based genotyping identified 58 AMR genes, including four blaCTX-M alleles (i.e., blaCTX-M-14, blaCTX-M-15, blaCTX-M-27, and blaCTX-M-55). Hierarchical clustering using selected phenotypic and genotypic characteristics suggested clonal dissemination of blaCTX-M genes. Out of 45 blaCTX-M gene-carrying strains, 7 harbored the gene in their chromosome. In one E. cloacae and three E. coli strains, chromosomal blaCTX-M-15 was located on insertions longer than 10 kb. These insertions were bounded by ISEcp1 at one end, exhibited a high degree of nucleotide sequence homology with previously reported plasmids, and carried multiple AMR genes that corresponded with phenotypic AMR profiles. CONCLUSION: Our study revealed the co-occurrence of ISEcp1-blaCTX-M-15 and multiple AMR genes on chromosomal insertions in E. cloacae and E. coli, suggesting that ISEcp1 may be responsible for the transposition of diverse AMR genes from plasmids to chromosomes. Stable retention of such insertions in chromosomes may facilitate the successful propagation of MDR clones among these Enterobacteriaceae species.

Distribution of amoebal endosymbiotic environmental chlamydia Neochlamydia S13 via amoebal cytokinesis.

We previously isolated a symbiotic environmental amoeba, harboring an environmental chlamydia, Neochlamydia S13. Interestingly, this bacterium failed to survive outside of host cells and was immediately digested inside other amoebae, indicating bacterial distribution via cytokinesis. This may provide a model for understanding organelle development and chlamydial pathogenesis and evolution; therefore, we assessed our hypothesis of Neochlamydia S13 distribution via cytokinesis by comparative analysis with other environmental Chlamydiae (Protochlamydia R18 and Parachlamydia Bn9 ). Dual staining with 4',6-diamidino-2-phenylindole and phalloidin revealed that the progeny of Neochlamydia S13 and Protochlamydia R18 existed in both daughter cells with a contractile ring on the verge of separation. However, in contrast to other environmental Chlamydiae, little Neochlamydia S13 16S ribosomal DNA was amplified from the culture supernatant. Interestingly, Neochlamydia S13 failed to infect aposymbiotic amoebae, indicating an intimate interaction with the host cells. Furthermore, its infectious rates in cultures expanded from a single amoeba were always maintained at 100%, indicating distribution via cytokinesis. We concluded that unlike other environmental Chlamydiae, Neochlamydia S13 has a unique ability to divide its progeny only via host amoebal cytokinesis. This may be a suitable model to elucidate the mechanism of cell organelle distribution and of chlamydial pathogenesis and evolution.

Corrigendum: Networking and Specificity-Changing DNA Methyltransferases in Helicobacter pylori.

[This corrects the article DOI: 10.3389/fmicb.2020.01628.].

Networking and Specificity-Changing DNA Methyltransferases in Helicobacter pylori.

Epigenetic DNA base methylation plays important roles in gene expression regulation. We here describe a gene expression regulation network consisting of many DNA methyltransferases each frequently changing its target sequence-specificity. Our object Helicobacter pylori, a bacterium responsible for most incidence of stomach cancer, carries a large and variable repertoire of sequence-specific DNA methyltransferases. By creating a dozen of single-gene knockout strains for the methyltransferases, we revealed that they form a network controlling methylome, transcriptome and adaptive phenotype sets. The methyltransferases interact with each other in a hierarchical way, sometimes regulated positively by one methyltransferase but negatively with another. Motility, oxidative stress tolerance and DNA damage repair are likewise regulated by multiple methyltransferases. Their regulation sometimes involves translation start and stop codons suggesting coupling of methylation, transcription and translation. The methyltransferases frequently change their sequence-specificity through gene conversion of their target recognition domain and switch their target sets to remodel the network. The emerging picture of a metamorphosing gene regulation network, or firework, consisting of epigenetic systems ever-changing their specificity in search for adaptation, provides a new paradigm in understanding global gene regulation and adaptive evolution.